Emergent Global Oscillations in Heterogeneous Excitable Media; 
The Example of Pancreatic (3 Cells 
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Using the standard van der Pol-FitzHugh-Nagumo excitable medium model I demonstrate a novel generic 
mechanism, diversity, that provokes the emergence of global oscillations from individually quiescent elements in 
heterogeneous excitable media. This mechanism may be operating in the mammalian pancreas, where excitable 
P cells, quiescent when isolated, are found to oscillate when coupled despite the absence of a pacemaker region. 
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I. EXCITABLE MEDIA 



An excitable element is defined by its response to a pertur- 
bation: whereas a small disturbance causes merely an equally 
small response, a perturbation above a certain threshold in 
amplitude excites a quiescent element that then decays back 
to quiescence during a refractory period in which it is un- 
responsive to further excitation. Such elements, when cou- 
pled to their neighbors into an assembly, become an ex- 
citable medium [[TJ . These have attracted an enormous amount 
of interest from different areas of science: the Belousov- 
Zhabotinsky reaction [^], plankton populations [0p, and the 
heart are just a few well-known examples. Oscillations 
in excitable media such as the heart are produced by forc- 
ing the medium from a pacemaker region; in other cases all 
individual elements of the medium oscillate when isolated: 
they are all pacemakers, and the medium is then oscillatory 
rather than excitable. However, this need not necessarily be 
so. In this paper I present a novel generic mechanism for the 
production of global oscillations; the introduction of diversity 
amongst the elements leads to the destabilization of the qui- 
escent state of an excitable medium and to the emergence of 
global oscillations even when each individual element of the 
medium is quiescent in isolation. I argue that an instance of 
such behavior may be found in a physiological example of an 
excitable medium without a pacemaker: the pancreatic (3 cells 
of mammals. 

In the mammalian pancreas are encountered structures 
where insulin is produced. These, the islets of Langerhans, 
consist of several thousand spherical cells clustered together 
and electrically connected via resistive gap junctions [ph. The 
vast majority of these cells, those that produce insulin in re- 
sponse to the level of glucose in the blood, are of a type known 
as (3 cells. In many mammals, including humans, the electrical 
potentials of (3 cells in an islet are found by experiment to cy- 
cle synchronously in slow oscillations termed bursts [Q]. On 
the other hand, an individual (3 cell when removed from the 
islet does not oscillate in this way [g], but rather is excitable. 
An islet of Langerhans, then, like the heart, is a physiological 
excitable medium. But in the pancreas, unlike in the heart, os- 
cillations of the excitable medium are not driven by pacemak- 
ing cells. The challenge is to understand the origin of these 
oscillations. Various detailed biophysical models of oscilla- 
tions in networks of pancreatic (3 cells have been proposed, in 



which the importance of the heterogeneity of individual cells 
U and the emergence of oscillatory behavior upon coupling 
nonoscillatory cells [10,11]] have been highlighted. Here I take 
a different approach in which I bring together these two ideas 
while simplifying the physics as far as possible to produce a 
minimal qualitative model for the phenomenon. I demonstrate 
that given diverse excitable elements, coupling these into a 
heterogeneous excitable medium can lead to the emergence 
of oscillations: that diversity is a generic mechanism for the 
emergence of global collective behavior not just in j3 cells, but 
in any heterogeneous excitable medium. 



II. VAN DER POL-FITZHUGH-NAGUMO MODEL 
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FIG. 1. An electronic excitable medium, (a) An element of the 
medium: the circuit and (b) the i-v characteristic of the nonlinear 
element nf. the cubic function i oc t> 3 /3 — v. (c) Resistive (diffusive) 
spatial coupling, illustrated for simplicity here in one dimension. 

To illustrate the mechanism, and having in mind its appli- 
cation to pancreatic (3 cells, I take an electronic circuit model 
as a caricature of a physiological excitable medium and as 
representative of excitable media in general. Each element 
of the medium is a circuit shown in Fig. [j](a). In physiolog- 
ical terms, the capacitor C represents a cellular membrane 
charged by E, characterizing ion pumps, and drained by a 
nonlinear resistance ni across C. This nonlinear element, a 
device with a range of negative resistance, could be a tunnel 
diode, for example, and should have the cubic i-v character- 
istic of Fig. |l](b). The inductance L models the finite switch- 
ing time of the ion channels in the membrane. The circuit is 
then mathematically described by the van der Pol-FitzHugh- 
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Nagumo equations [1121-115 



fj = — 7" (V? + v + f3rj). 



(la) 
(lb) 



Variables tji and rj are respectively proportional to the poten- 
tial difference across the nonlinear device ni and the current 
through the supply; parameter v is proportional to the poten- 
tial E, p to the resistance R, and 7 to the quotient L/C. The 
equations have an equilibrium point that is stable for |^| > S 
and unstable for \v\ < 3, where 
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so a circuit element is oscillatory for \v\ < 3 and excitable 
for \i>\ > 3 in the vicinity of = 5 [|l6|]. In physiological 
terms, oscillatory behavior corresponds to bursting < 3), 
and excitability to silent and continuously active cells; one 
case being v < — 3 and the other ^ > 3. 

Each of these elements is coupled to its nearest neighbors in 
one, two, or three dimensions to become an excitable medium. 
The coupling in biological and chemical excitable media is 
diffusive, though elastic excitable media that arise in electron- 
ics and rheology have also recently been considered [|l7 18). 
Either or both of Eqs. ([!]) may host a coupling term, depending 
on the medium being modeled. Here, to imitate a cellular ex- 
citable medium, the circuits are coupled resistively as shown 
in Fig. |l](c). This leads to a diffusive term in Eq. ( Jla] > 
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where the is represent the n neighboring elements, and k is 
the coupling strength or diffusion coefficient. In the contin- 
uum limit, the coupling term «X)i=iW'i ~~ VO becomes the 
Laplacian KV 2 ip; this is the classical van der Pol-FitzHugh- 
Nagumo model of an excitable medium, extensively analyzed 
in the excitable spiral-wave regime [|T^] . 



III. HETEROGENEOUS EXCITABLE MEDIA 

Consider now what happens if we introduce a diversity of 
parameter values for the different elements of the medium. 
The position of the stable equilibrium point for Eqs. ([!]) de- 
pends on the parameters v and (3, so if we introduce a spread 
of parameter values across the medium, we change the equi- 
libria of individual elements and the coupling between them 
K {ipi — VO w iH no longer be zero in the quiescent state. In 
such a heterogeneous medium, the dynamics of each element 
may be analyzed to a good approximation by treating the in- 
fluence of the rest of the medium as a type of external signal e. 
This is equivalent to setting e = k Y^i=i ~~ VO wnence the 
extra term e may be removed from Eq. (Bh by renormalizing 
v, and the element is modeled by Eqs. ([!]) with an effective v, 
v' = v — j3e; the medium is now oscillatory for \v — (3e\ < 5. 
Hence a small e can push an excitable element with \v\ > 3 



near to the excitable-oscillatory threshold over into the oscil- 
latory regime. This diversity mechanism is not specific to the 
van der Pol-FitzHugh-Nagumo model, but rather is generic; 
it can be applied to any heterogeneous excitable medium in 
which the position of equilibrium is parameter dependent and 
there is an oscillatory regime reachable in the extended pa- 
rameter space formed by an uncoupled element's parameters 
plus the extra term e. 
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FIG. 2. Dotted line: Temporal standard deviation of ip, a t 
(Eq. (Q)), shows the emergence of oscillatory behavior from a het- 
erogeneous excitable medium with increasing coupling n. Dashed 
line: Standard deviation of times t max of maxima of tp, a s (Eq. 
(^)), demonstrates increasing synchronization of the oscillations with 
increasing coupling /t. The numerical results are for a van der 
Pol-FitzHugh-Nagumo medium with /3 = 0.5, and 7 = 2, with 
4x4x4 elements randomly assigned v = 0.76 or v = —0.76, both 
of which are individually parameter values for which the medium is 
excitable rather than oscillatory, and thus quiescent unless excited. 

Let us examine the emergence of oscillations through di- 
versity in a simple example using the van der Pol-FitzHugh- 
Nagumo model. In Fig. |2|is plotted (dotted line) against cou- 
pling k the temporal standard deviation of ipj for a sample 
element j, 
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(4) 



for a heterogeneous van der Pol-FitzHugh-Nagumo medium 
consisting of elements randomly assigned the parameter val- 
ues v = 0.76 or v = —0.76, together with the other parame- 
ters j3 = 0.5, 7 = 2. Both of these parameter sets on their own 
produce excitable rather than oscillatory elements, which are 
quiescent without excitation, so at k = 0, a t is zero. But in 
Fig. H we see that as the coupling k between elements passes a 
threshold, cr t , which is measuring the temporal activity of the 
medium, increases from zero, meaning that the medium has 
spontaneously begun to oscillate; diversity has provoked the 
emergence of global oscillations from individually quiescent 
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elements. In this example, each element is connected in a cu- 
bic lattice with six neighbors. Elements throughout the lattice 
are randomly assigned one of two parameter values: v = 0.76 
or v = —0.76. On average, any element will find that half its 
neighbors share the same parameter value, and half have the 
other value. At quiescence, the coupling between those with 
the same parameter values is zero, while between those with 
different parameter values it is e = = 3kA?/>, 

since m = 3 is the average number of neighbors with differ- 
ent parameter values. In the van der Pol-FitzHugh-Nagumo 
model, the change in the equilibrium value of ip, Aip « — Av. 
For our example, v = ±0.76, which makes Av = 2v. In this 
case, then, the new effective v is 



|i/| = |i/-/3e| = |i/(l-6/3«)|. 



(5) 



We can see in Fig. g that the medium begins to oscillate when 
k w 0.06. For the numerical values (3 = 0.5, 7 = 2, 
\v\ = 0.76, used in Fig. g, |z/| = 0.62 when k = 0.06, which 
corresponds closely to the threshold for oscillation in a ho- 
mogeneous van der Pol-FitzHugh-Nagumo medium given by 
Eq. (Q) for these parameter values: \u\ = S = 31/96 y/7/2 = 
0.60412.... 

To simplify the above analysis as far as possible, I have 
considered a heterogeneous medium with just two states. The 
diversity mechanism works in the same way in a medium 
with a continuum of states, in which a proportion of elements 
— those whose uncoupled state falls within the parameter 
range for autonomous oscillation — may be intrinsic oscil- 
lators. However, such intrinsic oscillators do not act like the 
pacemaker region of a driven medium. Scattered randomly 
throughout a quiescent medium, they are neither sufficient nor 
necessary for global oscillations. Whether or not such intrin- 
sic oscillators are found in a heterogeneous medium depends 
on whether an oscillatory regime exists within the parameter 
range of heterogeneity for an element when the system cou- 
pling term e is zero. In the van der Pol-FitzHugh-Nagumo 
model this corresponds to whether the inequality \v\ < S is 
satisfied without coupling. For j3 cells, the majority of exper- 
iments have found no evidence for intrinsic oscillations, and 
that isolated cells are exclusively excitable: either silent, or 
continuously active rfjjj]. 



IV. SYNCHRONOUS AND ASYNCHRONOUS 
OSCILLATIONS 

The diversity mechanism applies very naturally to phys- 
iological excitable media, since homogeneous cells are a 
mathematical fiction. The underlying dynamics is such that 
there is a threshold for a quiescent heterogeneous medium to 
achieve criticality and cross the excitable-oscillatory bound- 
ary. This demands sufficient connectivity between neighbor- 
ing elements in terms of the number of connections (dimen- 
sionality of the system) and the coupling strength (diffusion 
coefficient), together with sufficient heterogeneity between el- 
ements. While the threshold-crossing mechanism of diver- 
sity is generic, what happens on the other side of the thresh- 
old once the medium has become oscillatory is not. In some 



oscillatory media, synchronous global oscillations are stable, 
whereas in others these are unstable to spatial perturbations 
leading to the formation of propagating fronts or other spatial 
phenomena 1 20 ] . The standard deviation across the medium 
of times t max (k) of maxima of ipk for elements k, 
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measures the spatial activity of the medium: the smaller 
this quantity, the greater the synchronization throughout the 
medium. In our example, a s plotted (dashed line) against cou- 
pling k in Fig. |] indicates that at the minimum coupling neces- 
sary for oscillatory behavior, the heterogeneous medium oscil- 
lates with little synchronization. As the coupling is increased, 
however, the synchronization rapidly improves, shown by the 
steep decay in the standard deviation of the maxima at larger 

K. 




FIG. 3. Plots of dispersion relations show de stabilization of syn- 
chronous global oscillations in the van der Pol-FitzHugh-Nagumo 
model for a range of v in the oscillatory regime. Here /? = 0.5, 
7 = 2, k = 0.09: shown are curves for v — 0.592 (dotted), 
v = 0.596 (dashed), v = 0.600 (solid), and v = 0.604 (dot-dashed). 
The maximal Floquet exponent A gives the growth rate per period of 
a perturbation of wavenumber q WM- For v = 0.596 and v — 0.600, 
the dispersion relation \ q as a function of q is positive for a range of 
q, meaning that there is a set of wavelengths for which perturbations 
grow exponentially. Outside this range of v values, the dispersion 
relation is never positive, implying linear stability of synchronous 
global oscillations under perturbations of all wavelengths. 

If the synchronized state is an attractor for the homoge- 
neous medium, then we might expect the behavior in the het- 
erogeneous case to reflect this. Let us consider for a moment 
the related area of synchronization, phase- and frequency- 
locking, or entrainment of coupled oscillators, a vast field 
of study initiated by Huygens with his observations of syn- 
chronization of two pendulum clocks coupled by a common 
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mounting [21- 23 1 . Winfree p4| ] showed that synchroniza- 
tion emerges in a population of heterogeneous oscillators as 
coupling exceeds a critical threshold, in a manner reminiscent 
of a thermodynamic phase transition. Following this work 
Kuramoto pa-El] developed his theoretical model whose 
tractability helped to advance significantly the study of cou- 
pled heterogeneous oscillators. This undergoes two transi- 
tions as the spread of natural frequencies is reduced, or the 
coupling is increased: first comes the onset of partial syn- 
chronization, which is followed by complete synchronization 
even with some residual heterogeneity among the elements. 
More recently a physically realizable version of the Kuramoto 
model has been proposed: an array of heterogeneous Joseph- 
son junctions Pq], Strong coupling between oscillators al- 
lows modification of the amplitude as well as the phase of an 
oscillator, which can give rise to novel phenomena such as 
amplitude death p^-pll] (here however we are in the weak 
coupling regime). Finally, if synchronous global oscillations 
are unstable in the homogeneous case, the addition of a cer- 
tain amount of heterogeneity can even lead to an increase in 
synchronization, as studies with locally coupled limit-cycle 
oscillators have found [32,331]. 
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FIG. 4. Dotted line: Temporal standard deviation of ipj, a t (Eq. 
(^)), illustrates the transition to oscillatory behavior in a homoge- 
neous excitable medium with decreasing v. Dashed line: Standard 
deviation of times t max of maxima of ip, a s (Eq. (^)), shows the 
destabilization of synchronous global oscillations in a range of v 
just above the excitable-oscillatory threshold. The inset highlights 
the narrow range 0.588 < v < 0.608 in which the destabilization 
occurs. Numerical results are for a van der Pol-FitzHugh-Nagumo 
medium with 4x4x4 elements with (3 = 0.5, 7 = 2, and k = 0.09. 
From Eq. (^), the medium is oscillatory for \v\ < 3 = 0.060412 .... 

We can obtain a theoretical understanding of synchroniza- 
tion in a homogeneous oscillatory medium by calculating an- 
alytically the Floquet exponents that indicate the linear stabil- 
ity of limit cycles against spatial perturbations. This I have 
done for the van der Pol-FitzHugh-Nagumo model in Fig. [| 
where I show a family of dispersion relations for different val- 



ues of v. A positive value for the maximal Floquet exponent 
illustrates the destabilization of synchronous global oscilla- 
tions of the medium for a range of v for a particular value 
of j3. In the case shown, for which, as in Fig. 0, /3 = 0.5, 
7 = 2, the dispersion relation X q is positive for a range of 
wavenumbers < q < q c for v values close to the excitable- 
oscillatory threshold of the medium. For these v values then 
synchronous global oscillations are unstable to perturbations 
of wavelengths 1 jq c < A < oo when these wavelengths fit 
within the system size. For all other v values in the oscillatory 
regime global synchronous oscillations of the homogeneous 
medium are linearly stable. This is demonstrated numerically 
in Fig. Q in which are plotted against v the measures of ac- 
tivity and spatial synchronization at and a s for a homoge- 
neous van der Pol-FitzHugh-Nagumo medium with (3 = 0.5, 
7 = 2, and k = 0.09. The dotted line er t shows that the 
medium is oscillatory up to v = 0.604, then excitable, as Eq. 
^) demands. The dashed line a s is positive near v = 0.6, 
showing that global oscillations are unstable for a range of v 
values; outside this range the standard deviation a s is zero, as 
the medium oscillates synchronously. The inset highlights the 
range of v in which destabilization occurs. 

In a homogeneous van der Pol-FitzHugh-Nagumo ex- 
citable medium we have seen that there are parameter ranges 
in which synchronous global oscillations are stable, and oth- 
ers in which they are not. The behavior of the homogeneous 
case helps us now to understand the heterogeneous medium. 
Figure || displays poor synchronization of the medium imme- 
diately following the emergence of oscillations, which then 
improves with increasing k: a s peaks, then declines rapidly 
before leveling off. This can be seen as a combination of 
the intrinsic dynamics of the van der Pol-FitzHugh-Nagumo 
medium together with the synchronizing effect of coupling. 
Upon the emergence of oscillations in the heterogeneous 
medium, the effective v, v' (Eq. (Q)), is within the range of 
values near to the excitable-oscillatory transition for which 
synchronous oscillations are unstable. As the coupling k in- 
creases, v' drops below this range and the synchronization im- 
mediately improves. The further slower decrease in a s is due 
to the additional synchronizing effect of coupling. 

What might be the relevance of the synchronization of 
global oscillations for the cells of the pancreas? The pan- 
creas may be contrasted with another physiological excitable 
medium: the heart. While in the heart synchronous oscil- 
lations are vital to the survival of the organism — the un- 
synchronized state of the fibrillating heart is fatal if not im- 
mediately resynchronized with an electric shock — it is not 
obvious why this should be so for j3 cells. Although in hu- 
mans and in mice the oscillations are synchronous, in other 
mammals there is less evidence for this. While physiologists 
may yet provide a biological rationale for synchronization, it 
may be that it is not a physiological necessity but simply a 
byproduct of the emergence of oscillations: the pancreatic [3 
cells of some species may be operating in a parameter range 
in which synchronous global oscillations are stable, while in 
others, less investigated up to now, spatial patterns may be 
found in the oscillations. 
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V. CONCLUSIONS 

I have argued that diversity acts at a fundamental level in the 
dynamics of heterogeneous excitable media to produce global 
oscillations. Another related theme of study has subjected ex- 
citable systems to various types of forcing. A periodically 
forced van der Pol-FitzHugh-Nagumo element shows behav- 
ior similar to a driven oscillator: phase locking, quasiperiod- 
icity, period doubling, and chaos |53]. A quiescent excitable 
element may be excited by driving with a combination of a pe- 
riodic subthreshold signal plus noise [ ]35| ] or with an aperiodic 
subthreshold signal plus noise [f36[], phenomena which have 
been termed stochastic resonance, or with noise alone J37|], 
when the phenomenon has been termed coherence resonance. 
Here, we have seen that even without any external forcing, ei- 
ther periodic or stochastic, a heterogeneous excitable medium 
can become self-excited to produce global oscillations. 

To what extent is this general analysis applicable to pan- 
creatic (3 cells? (3 cells are diffusively coupled excitable 
elements, although more complicated than van der Pol- 
FitzHugh-Nagumo elements, with more internal variables and 
parameters. Their heterogeneity has been made manifest in 
studies demonstrating differing rates among cells of insulin 
synthesis and secretion [|38|-|42[|, of glucose metabolism [p3|], 
of changes in calcium concentration [ p4p5| ], and of changes 
in electrical activity [Q. Physiologists have suspected the 
importance of this variability [^7|]. (3 cell oscillations, or 
bursts, are more complex than the simple oscillations of a van 
der Pol-FitzHugh-Nagumo type relaxation oscillator, but the 
basic mechanism that diversity provokes the emergence of os- 
cillations remains the same, and provides a qualitative expla- 
nation for the emergence of global oscillations in the (3 cells 
of the mammalian pancreas. Heterogeneity is the norm in bi- 
ological excitable media, so there may well be other instances 
awaiting discovery of this mechanism in operation. 
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